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ABSTRACT 

This article presents the application of a state feedback, Pole Placement temperature tracking 
control formulation leading to a comprehensive Single input, multiple output (S1MO) structure for 
cooling of a steel plate by an impinging air jet. . Following a semi discrete control volume formulation 
of the 1-D transient heal conduction: the state space model of nodal variation of metal temperature 
with lime has been arrived. Stable temperature track with respect to the reference temperature track 
was obtained by using ilie state feedback control algorithm. To solve this prescribed state feedback 
problem, an efficient pole placement technique was implemented in two ways i.e.. initially, the fixed 
poles and then variable poles in the complex plane. The control m a t i imula d in the M ' 

and SIMULINK environments. The state feedback pole placement technique was found efficient in 
controlling the parameters in the given application. It was observed that the response of the non-linear 
system is sensitive to linearization lime interval. Better control is implemented by increasing the 
frequency of adjustment of the closed-loop poles. 

KEYWORDS: Feedback Control, Jet Impingement Cooling, Pole Placement Technique, Slate Space 
Method, Single Input Multiple Output Problem. Fixed Poles, Variable Poles, Temperature Control 



Nomenclature: 

K: The feedback gain matrix t: Time (s) 

A: The coefficient matrix in the state space equation 
Ai. The coefficient matrix in the state space representation for i th model A : 

The mean coefficient matrix B: The coefficient matrix of control 

input T] : Temperature at the top surface of the steel plate (K) 

B ; : The coefficient mattix of control input or i th model 
X: The matrix of state variables X : The 

matrix of mean state variables X : The matrix of state rates 

a e : The perturbation in strain rates 
A :The eigen values X ' : The error matrix [X — X J 

T, : The demand (K/s) P, (S): The 

ith number of transformable models a: The control input (Strain rate) 

( S ~ ) a : The mean strain rate (mean control output) 

A-BK: The feedback matrix T: The 

transformation matrix I : The identity matrix 

Z: The assumed state variable matrix r: The 

number of operating conditions 

I. INTRODUCTION 

Heat treatment involves the method of controlled heating, soaking and cooling temperature cycles that 
are applied to a material or substrate, to impart superior metallurgical properties to it. Heat treatment has wide 
application and is most abundantly used in the field of glass, ceramic, iron and steel production industries. 



Model bused numerical stale feedback... 

The dynamics of heat transfer systems are highly nonlinear. Irrespective of nonlinearity in of heat 
transfer systems, air jet impingement cooling systems also have model uncertainties which includes large 
changes in the air flow rate and variation of parameters, the external disturbances, leakages, and friction. To 
overcome this problem in order to gain profits, process designers are creating designs and modelling in complex 
non-linearity regions where process controllers can face stiff the challenges. In model predictive control, the 
controller action is the solution to a constrained optimization problem that is solved numerically on-line [1,2, 
3]. In contrast, differential geometric control is a direct synthesis approach in which the controller is derived by 
requesting a desired closed loop response in the absence of input constraints [4, 5]. In Lyapunov -based control, 
closed loop stability plays a central basic role in a controller design. [6, 7, 8]. Kravaris and Daoutidis [9] 
presented a non-linear state feedback controller for second order non-linear systems. Niemiec and Kravaris [10] 
proposed a systematic procedure for the construction of statically equivalent outputs with prescribed 
transmission zeros. Kanter et al. [11] developed nonlinear control laws for input constrained, multiple -input, 
multiple-output, stable processes. They addressed the nonlinear control of the processes by exploiting the 
connections between model- predictive control and input-output linearization. Kravaris et al. [12] presented a 
systematic method arbitrarily assigning the zero dynamics of a nonlinear system by constructing the requisite 
synthetic output maps. The method requires solving a system of first-order, nonlinear, singular PDEs. Mickle et 
al. [13] developed a tracking controller for unstable, nonlinear processes by using trajectory linearization. 
Tomlin and Sastry [14] derived tracking control laws for non -linear systems with both fast and slow, possibly 
unstable, zero dynamics. Vander Scaft [15] developed a nonlinear state feedback H ™ optimal controller. 
Devasia et al. [16] and Devasia [17] introduced an inversion procedure for nonlinear systems that constructs a 
bounded input trajectory in the pre-image of a desired output trajectory. Hunt and Meyer [18] showed that under 
appropriate assumptions the bounded solution of the partial differential equation of Isidori and Byrnes [19] for 
each trajectory of an ecosystem must be given by an integral representation formula of Devasia et al. [17] Chen 
and Paden [20] studied the stable inversion of nonlinear systems. 

The present work uses the same continuous-time, model predictive control framework. Conceptually as 
the amount of nonlinearity is very less so the model does not require any linearization of the process equations. 
Here the nonlinear terms are considered as the strain rate in the convective heat transfer correlation. The 
controller introduced here is obtained by requesting desired responses for the state variables. The nonlinear state 
feedback is derived by minimizing a function norm of the deviations of the controlled outputs from linear 
reference trajectories with orders equal to the (output) relative orders. In the part of the work described in this 
chapter, we come up through a systematic evaluation with a process model that provides the starting point of 
evolving a model-based control. At the same time, we analyze the performance of a model-based nonlinear 
control with the applications of an efficient pole-placement technique. The focal point of the present analysis is 
a detailed modeling of the heat transfer and applications of efficient techniques (pole -placement) to control the 
state variables influencing the process. Moreover this chapter presents a nonlinear control method that is 
applicable to stable and unstable processes. The closed loop stability is ensured by forcing all the process state 
variables to follow their corresponding reference trajectories. The proposed control system includes a nonlinear 
state feedback and pole placement technique in both the forms of fixed poles and variable poles. 

1.1. Pole placement techniques 

"Pole placement" is a conventional design methodology for linear time-invariant control systems. The 
method is based on the fact that several performance specifications can be met by using dynamic output 
feedback to adequately placed closed loop poles in the complex plane. . An extension of the classical pole 
placement problem is the regional pole -placement problem, in which the objective is to place closed-loop poles 
in a suitable region of the complex plane. In many real cases, the model uncertainty reflects on the parameters 
of the plant, which has motivated extensive research efforts in parametric robust theory [21, 22, 23]. Pole 
placement methods can be divided mainly into two categories such as output feedback and state feedback 
methods. Usually dominant pole design techniques are employed as it is not always possible to assign all the 
poles of a system and thus one has to relax the design requirements and consider a partial pole placement that 
will give satisfactory overall performance. The pole placement problem is much easier to handle when the 
design is done in the state space domain and many techniques exist for this particular case. In this work, the 
plant is represented by proper state-space equation. Robust pole placement problem is formulated as the 
problem of robustly assigned poles in a region in two different ways i.e., by assuming some appropriate fixed 
poles and variable poles in terms of the Eigen values of the ( A — BK ) matrix. Where 'K' is the gain. Similar 
conditions hold in [24] and [25]. 
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The main focus of this work was to present a more efficient way of approaching Ackerman's formula 
to find the gain matrix so as to achieve the Eigen poles. Apart from this, the technique presented here was 
utilized towards significant progress for the case of nonlinear Single-input and multiple-output (SIMO) systems. 
. As only the temperatures of the plate are controlled, typically the plate may remain insignificantly cool or the 
desired cooling temperature trajectory is followed backwards. That is why we focus on state tracking instead of 
output tracking. To facilitate the design of controller, we first investigate the pole placement problem for SIMO 
systems using dynamic gains, when all the state variables are available for feedback. Both the local and global 
pole placement control has been accomplished in the MATLAB and SIMULINK environments. In both the 
cases the results are shown in the comparative mode between the demand and controlled states. 



II. THE DESIGN PROBLEM 

In this work, we consider the modeling to evolve a dynamic model-based control of the air jet 
impingement cooling of steel plate. Here we propose the design procedure based on state space pole placement 
techniques for the state feedback control design of jet impingement cooling system. The design of State-Space 
models is not different from that of transfer functions in that the differential equations describing the system 
dynamics are written first. For state-Space models, models, instead, the equations are arranged into a set of first 
order differential equations in terms of selected state variables, and the outputs are expressed these same state 
variables, Because the elimination of variables between equations is not an inherent part of this process, state 
models can be easier to obtain. The State models are directly derived from the original system equations. The 
standard form of the State-Space model equation is: 

X(t) = AX(t)+ Ba(t) (State Equation) 
Here 'x' is the State Vector, the vector of state variables, a is the control variable and 'A' is the system 
coefficient matrix. The matrices A and B are controllable. A fixed dynamic controller is designed i.e., 

a = B ^[AX-X] 

The pole placement design of 'B' in the form of full state feedback becomes equivalent to design 'K' 
where, 

a e = -KX e (Controller or Control Law) 

(1) 

For stability all the Eigen values of 'A' must be such that the real part is negative. Let us assume, 
Z = TX 

(2) 



The equation becomes, 
TX = TAX + TBa 
Z = TAT ^Z + TBa 



,- X = T _1 Z) 



Z = AZ + Ba ; 

Where A = TAT -1 and B = TB 

Let us consider the Eigen values of equation (3) be X , 
Then, | A - Xl| = 

|a-3j| = o=>;i=X 
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Figure 1: State Feedback Configuration 
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Since, X = X , the transformation "I" does not alter the system stability. The above diagram represents the 
corresponding control structure that evolves in the form of a feedback structure. In contrast to traditional 
feedback structure, K varies with time. Then feeding this input a(t) back into the system, we obtained, 

X(t) = (A-BK)x(t) 
(4) 

III. FULL STATE FEEDBACK DESIGN 
There are several motivations for formulating the design procedure in state space. One motivation is a 
clear exposition of the controller structure. Another is that the controller structure points possibilities of 
reducing the order of the compensators using state space technique [26]. A third motivation is that the state 
space formulation may accommodate design procedures for systems in which the number of inputs or outputs is 
less than the number of operating conditions. 

Let us consider the system P has r > 1 operating conditions with models P ; (S), P 2 (S), P r (S). The state 

space representations of the models P; (S), i = 1,2 . . . r, are: 

X = A i X + B i a 

(5) 

Where the state X e R n and the inputa e R m . The pairs [A. i ,Bj J , i = 1,2. . .r are controllable. Let us design 
a fixed static controller a = — KX such that the poles of the closed loop system are ^.(Aj — BjKJ =A i , i = 
1 . . .r, where <T i are the desired Eigen values for the i lh operating conditions. For full state to have a solution, it 
is necessary that the pairs (A ; , -B, ), i= 1 . . .r, be controllable. Closely related to full state problem is the design 
of a fixed gain matrix K to stabilize all the operating conditions. The solution to problem exists a K such that: 

det(sl - (A. - BjK)) = Pi (s),i = l,...,r. 

(6) 

The nonlinear pole-placement equation (6) can be considered as Bass-Gura pole-placement formula [27], 
extended to multiple operating conditions. 

IV. NUMERICAL METHODS FOR THE SINGLE POLE PLACEMENT PROBLEM 

There are several methods for the pole placement problem some of the well known theoretical 
formulae, such as the Ackermann's formula, Bass-Gura formula etc.. The primary reason is that, these methods 
are based on transformation of the controllable pair (A ( , Bj ) to the controller companion form. 
The computationally viable methods for pole placement are based on transformation of the pair (A^Bj J to 
the controller Heisenberg form or the matrix A to the real Schur form (RSF), which can be achieved using 
orthogonal transformations. The pole placement design in the present study is accomplished in two ways, that is 
i) Assigning the fixed poles arbitrarily in the complex plane, 

ii) Assigning the variable poles in the complex plane. 
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Figure 2: Block Diagram Showing Pole-Placement Control Model for Jet Impingement 
Cooling of Steel Plate (Temperature tracking Control). 
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V. STATE FEEDBACK CONTROL BY FIXED AND VARIABLE POLE 

PLACEMENT DESIGN 

As a first attempt, a non-linear model-referenced state-feedback tracking control by fixed (local) pole 
placement for a jet impingement heat transfer system was implemented. The closed-loop poles were adjusted 
arbitrarily at each of the time intervals based on the open -loop poles of the system. The arbitrary pole 
assignment was performed in the SIMULINK environment. The poles were selected on the basis of trial and 
error method, for different simulation times. For each time segment the controlled output response and the 
reference output track were identified and compared. The results obtained in the SIMULINK environment for 
different time segments and time steps adopting the fixed pole placement control of the temperature track are 
discussed in detail in the later sections. After picking the feedback gain matrix K to get the dynamics to have 
some nice properties (i.e., Stabilize A) X i (A) — > 'k i (A — BK) , the selected poles were placed in the complex 
plane. Recalling that the characteristic polynomial for this closed-loop system is the determinant 
of(sI — (A — BKJJ. Since the matrices A and B*K are both 3 by 3 matrices, there were 3 poles for the 
system. By using full state feedback we placed the poles at different places arbitrarily. The control matrix K, 
which gave the desired poles, was found by using the MATLAB function PLACE. For SIMO systems, PLACE 
use the extra degrees of freedom to find a robust solution for K i.e., minimizes sensitivity of closed loop poles 
to perturbations in A and B. 

Then, the pole placement design procedures have been executed for the state feedback design of systems with 
single operating condition. To facilitate the design, poles varying with time were introduced. The same 
controller using strictly the state space technique was used for simulation i.e., a = — KX . This approach 
provided a clear exposition of the controller structures. The closed-loop poles were adjusted at each of the time 
intervals based on the open-loop poles of the system. The feedback gain was developed based on the plant 
models and the corresponding set of closed-loop poles. The response of the non-linear plant to this controller 
was presented for the full state feedback. In order to target a steady state of the system the Eigen values of (A- 
BK) should be such that their real parts are negative. They were chosen in the vicinity of the open-loop poles of 
the system. If any of the open-loop poles have positive real part, the real part of the corresponding closed-loop 
poles were chosen to be zero. The time variant Eigen values of the matrix (A-BK) were computed then 
computed. 

VI. RESULTS AND DISCUSSIONS 

The main focus of this work was to evolve a model based closed loop control analysis with the 
application of an efficient pole placement technique. On the way to establish the necessary control model, we 
captured an important jet impingement cooling of a steel plate. The simulations were carried out in the 
MATLAB and SIMULINK environment. Different sets of results were generated for different set of poles. The 
results discussed below are both for fixed closed loop pole placement and variable pole placement problems. In 
the closed loop control phenomena, the important task was evaluation of elements of the feedback gain matrix 
(K). 

6.1 Closed- loop control results for fixed pole placement technique 

is observed that with the application of fixed poles such as two complex conjugate poles and one real 
pole, where all the real poles are negative i.e., poles are considered more in the left half of the complex plane, 
the cooling rate controls upto -lK/s and the nature of controlled temperature track is almost coincides over the 
reference temperature track as shown in the figures 3 and 4. The poles selected initially are —0.2 + 1/ 
(Complex Conjugate) and -100 (Real). Several trials were carried out by changing the pole locations and plots 
were generated for the temperature track on the top surface of the steel plate to realize the controlled cooling 
rate. 
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Figure 4: Effect of linearization time interval on tracking control of Plate temperature by 

assigning a fixed set of poles at Time Steps of 100 and Time segment of 1 second. 
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6.2 Closed- Loop control results for variable pole placement technique 

This technique applies the Eigen values of the matrix (A) where, the matrix A is a 3 by 3 matrix, so 
there will be 3 poles for the system. Using the MATLAB function eig, the Eigen values of the (A) matrix are 
evaluated to evolve the feedback gain matrix (K). The problem was solved as a closed loop feedback control 
model using time dependent poles. For the positive values of the poles, the stability of the system was found 
marginal. Thus, at each linearized time segment the appropriate poles were assigned and the corresponding gain 
matrix was computed. The total computation has been performed by the help of SIMULINK tools. The effects 
of linearization time segment on tracking control of Plate temperature by assigning variable poles at different 
time locations are presented in the figures 5, 6, and 7 respectively. The only peculiarity found in this case was 
that, the response was oscillatory in nature. These results clearly predict that, one of the dominant poles tend to 
more negative side, while other two become less effective. 



- 


-- 






° e 


r--- 












1 = 


~ 


. - : - 


-i < = 


1 












" 




" ?isVs>Hl" 






(A-EK> 











Figure 5: Effect of Linearization Time Interval on Tracking Control of Plate Temperature by Assigning a 
Set of Variable Poles at Different Time Locations at Time Steps of 100 and Time Segment of 0.1 second. 




Figure 6: Effect of Linearization Time Interval on Tracking Control of Plate Temperature by Assigning a 
Set of Variable Poles at Different Time Locations at Time Steps of 1000 and Time Segment of 1 second. 
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Figure 7: Effect of Linearization Time Interval on Tracking Control of Plate Temperature by A 
Set of Variable Poles at Different Time Locations at Time Steps of 100 and Time Segment of 1 second. 
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VII. CONCLUSION 

A model based simulation and control of the temperature track in jet impingement cooling process of a 
steel plate has been developed that incorporates the implementation of an efficient pole placement technique in 
two ways i.e., one by the local poles and other by the global poles. The nonlinear control model developed here 
can be used to evaluate the performance in terms of physical properties and metallurgical aspects. The 
conclusions being obtained from the work are: 

1. Pole placement design procedures have been proposed for the state feedback design of systems with single 
input condition. To facilitate the design, the concepts o coefficient matrices (A and B) and local poles are 
introduced. 

2. It is seen that the response of the non-linear system is sensitive to linearization time interval. Better control 
is implemented by increasing the frequency of adjustment of the closed-loop poles. 

3. In the present algorithm, the mean solution is utilized at discrete intervals (Piece-wise linearized time 
intervals) of time for evaluation of the coefficients. 

It is observed from both the fixed pole placement and variable pole placement that, in case of fixed 
pole placement, the simulation time is limited uptol.4 seconds. Within this simulation time, the demand 
temperature track and the response are coinciding to each other with very negligible error. While, in case of 
variable pole placement, the simulation time is large for the same linearized time intervals and percentage of 
error is significant compared to the first case. 
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